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O '■ Abstract 



Lattice QCD should allow quantitative predictions for the heavy quark 
physics from first principles. Up to now, however, most approaches have based 



Oh! 

,^ . on the nonrelativistic effective theory, with which the continuum limit can not 



be taken in principle. In this paper we investigate feasibility of relativistic ap- 
proaches to the heavy quark physics in lattice QCD. We first examine validity 
of the idea that the use of the anisotropic lattice could be advantageous to 
control the mqa corrections. Our perturbative calculation, however, reveals 
that this is not true. We instead propose a new relativistic approach to handle 
heavy quarks on the isotropic lattice. We explain how power corrections of 
mqa can be avoided and remaining uncertainties are reduced to be of order 

(oAqcd)^. 
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I. INTRODUCTION 

Weak matrix elements associated with B mesons are essential ingredients to determine 
Cabibbo-Kobayashi-Maskawa matrix elements. In principle lattice QCD provides the oppor- 
tunity of evaluating these matrix elements from first principles. However it is still difficult 
to simulate the b quark with high precision on the lattice. The main source of systematic 
errors originates from the fact that the b quark mass in the lattice unit is large: rriba ~ 1 — 2 
in the quenched approximation and rriba ~ 2 — 3 in full QCD with current accessible compu- 
tational resources. In order to control large rriQa errors, several ways have been proposed so 
far: A static approximation with rriQ — i> oo [|ll, a nonrelativistic QCD 0, a nonrelativistic 
interpretation applied to results with the Wilson/Sheikholeslami-Wholert (SW) action |^ 
and an anisotropic lattice with finer temporal lattice spacing at while keeping the spatial 
one ttg modest @] . Although the b quark can be directly simulated with any of the last three 
approaches, only the last one has the advantage that we can take the continuum limit, which 
is a fascinating feature stimulating our interest. 

Practical effectiveness of the anisotropic lattice is transparent: with finer temporal lattice 
spacing time evolutions of all kinds of correlation functions become milder, which benefits us 
the better signal-to-noise ratio. On the theoretical side, our interest exists in the use of the 
anisotropic lattice to control mqa errors. If they are restricted to only powers of mqat, they 
can be made small by the anisotropic lattice with smaller at. Indeed mqas corrections can be 
removed at the tree level |]5|. Our main concern, however, is whether mqas corrections could 
revive perturbatively or nonperturbatively. Up to now no one has successfully eliminated 
the possibility that mqas corrections can appear beyond the tree level. Without the proof 
given the anisotropic lattice is no better than the isotropic one, where one can also eliminate 
mqa corrections at the tree level, in terms of controlling m,qa corrections. 

The first part of this paper is devoted to a one-loop calculation of the quark self energy 
on the anisotropic lattice. We analyze its 0{g^a) terms to examine the possibility of ap- 
pearance of 0{g'^m,qas) contributions. Our results strongly suggest that one- loop radiative 



corrections allow the revival of mqas corrections. Since we find that the anisotropic lattice 
is not theoretically advantageous any more, in the second part of this paper we propose 
a new relativistic way to deal with the heavy quarks on the isotropic lattice, analyzing 
rriQa corrections carefully. We discuss cutoff effects of the heavy quark system following 
the on-shell improvement programme |§|,0|. An important finding is that leading cutoff ef- 
fects of order (mQa)"' can be absorbed in the definition of renormalization factors for the 
quark wave function and mass, so that the leading rriQa correction is reduced to be of order 
{rriQa)'^ a Aqcb- After removing remaining leading cutoff effects of OllmQa)"^ a Aqqd) with 
parameters in the quark action properly adjusted, we are left with only 0((aAQCD)^) errors, 
which are expected to be fairly small. We also provide a nonperturbative method to control 
the 0((mga)") corrections involved in renormalization factors. 

This paper is organized as follows. In Sec. |I| we present a one-loop calculation of the 
quark self energy on the anisotropic lattice and discuss the possibihty of the revival of rriQas 
corrections. In Sec. ^ we propose a new relativistic approach to handle the heavy quarks 
on the isotropic lattice avoiding large rriQa corrections. Our conclusions are summarized in 
Sec. 0. 

II. ANISOTROPIC LATTICE 

A. On-shell improvement on the anisotropic lattice 

In order to obtain a generic form of the quark action allowed on the anisotropic lattice, 
let us make the operator analysis according to the Symanzik's improvement programme 
PJ^. The lattice theory is described by a local effective theory as 

'S'eff = 5*0 + ^ a^ j d'^xc4+k,i{9)04+k,i{^), (1) 

k>l,i •' 

where 5*0 denotes the continuum action. Ok,i{x) is a local composite operator with k dimen- 
sions, which consist of quark mass, quark fields and link variables. These higher dimensional 



operators must respect symmetries on the lattice such as the invariance under gauge, par- 
ity and charge-conjugation transformations and discrete rotations. The coefficient Ck^i is a 
function of the gauge couphng g to be determined perturbatively or nonperturbatively. 

Symanzik's improvement programme was originally designed to reduce cutoff effects order 
by order for on-shell and off-shell Green functions. However, we are interested in only on-shell 
quantities such as hadron masses and matrix elements which require correlation functions at 
non-zero physical distances. Here it would be better to consider the on-shell improvement 
procedure that is much simpler but restricted to on-shell quantities [§. 

Under the requirement of various symmetries on the lattice, we find the following set of 
operators with dimension up to five: 

dim.3 -.O'^ix) = q{x)q(x), (2) 

dim.4 :04„(x) = q{x)-foDoq{x), (3) 

0',,ix) = Y.q{xhiD,q{x), (4) 

i 

dim.5 -.O'.^ix) = q{x)Dlq{x), (5) 

o',,{x) = Y.qi^)DH^), (6) 

i 

<^5cix) = i H g(x)aoiFoig(x), (7) 

i 

O'^di.^) = i Y. q{x)aijFijq{x), (8) 

O'.M = J2^i^)boDo,liD^]q{x), (9) 

i 

where D^ is the covariant derivative on the lattice and afj,^ and F^i^ are defined as a^^ = 
[7^,7y]/2 and igF^^, = [D^,D,^]. The subscript denotes the time component, while i,j = 
1,2,3 space components. In terms of these operators, a general form of the quark action on 
the anisotropic lattice is given by 



ST = E 



40^3:) + E 4c?L(^) + E 4A 



5i 



X 



(10) 



where C3, . . . , Cgg are functions of the bare gauge coupling g, the bare quark mass niQ and 
the time and space lattice spacings at, a^. Since degrees of freedom of C3 and €4^ can be 



absorbed in the renorinalization of the quark mass and the wave function respectively, we 
choose C3 = mo and c'^^ = 1 for convenience. We also find O'^^ and Og^ are related to 
other operators by using the classical field equation, and hence Cg^ and Cgg can be set by 
hand. We eliminate Cgg by choosing Cgg = and employ O'^^ to avoid species doubling with 
Cg^ finite. After all remaining parameters are 04^, Cg^, Cg^ and Cg^, which should be tuned 
(nonperturbatively) in order to remove 0{at^s) discretization errors, so that we are left with 
only discretization ambiguities of 0{a^J. This point should be stressed, since previous 
papers P,|] claim that only three parameters 04^, Cg^ and Cg^ are enough to be tuned for the 
0{at^s) improvement. 

B. Quark and gauge actions 



According to the discussion in the above subsection a general form of the quark action 
on the anisotropic lattice is given by: 



Sq = atal Y^ g(5 



at 



^oDo + iyJ2 7. A + mo - ^r (Do' + r/ E A' 



tg 



-at—r <CEil+r])J2 ^oi.Foi + cbV J2 ^ij^i^ 

i ij 



q{x) 



where parameters in eq. (|T0|) are rewritten as 



C46 = ^, 



1 


r 


4a = 


"2' 


4 = 


rrj 
2' 


4 = 


-|rc^(l + r/) 


4 = 


9 

--rcBV- 



(11) 



(12) 
(13) 
(14) 
(15) 
(16) 

Here the Wilson parameter r is taken arbitrary by hand as mentioned in the previous 
subsection. Note that the quark fields q, the covariant derivatives Dq and Di, the bare quark 
mass mo, the gauge field strengths Foi and Fij are dimensionful. They are transformed into 
dimensionless quantities by 



q = al/^q, 



Dq = atD. 



0; 



Di = ttsDi, 



niQ = atmo, 
Fqi = atagFoi, 



F- = n'^F- 



(17) 
(18) 
(19) 
(20) 
(21) 
(22) 



With this transformation, eq. ([TlD is rewritten as 



Sq{x) =J2^(^) 



^ i ^ i 



ig \ (1 + 7]) 






U 



q{x) 



with ^ = tts/at the anisotropy parameter. 

From eq. (|Tl|) we find the inverse free quark propagator 



(23) 



z/ 



(^tSg ^{p) = z7osin(poai) + -iY^%sm{pias) + moat 

^ i 

+r(l - cos(poat)) + 7^ Yli^ ^ cos(pia^)), 

? i 

and relevant vertices for the present calculation 
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-CEgT (Tjosm (poctt - goflt) cos 



-csgT — E cTijSin (pj-a^ - g^a^) cos 



(24) 



(25) 
(26) 
(27) 
(28) 
(29) 

(30) 



where p is for the incoming momentum into the vertex and q for the outgoing momentum. 
T"^ {A = 1,...,N^ -1) is a generator of color SU(A^c)- 

For the gauge part we take the standard Wilson action on the anisotropic lattice: 



c _ 2iV 



^j:(^-7^T^^iUUn) + UUn)) 



2N, 



7T.(^-7^T^<U,An) + U},{n)) 






(31) 



with 



U,,{n) = U,{n)U,{n + il)Ul{n + i>)Ul{n). (32) 

The gluon propagator is given by 

in the Feynman gauge with the fictitious mass A^, which is introduced to work as the infrared 
cutoff. 

C. Tree-level analysis on the quark propagator 

Expanding the inverse free quark propagator in eq. (0) up to 0{at) we obtain 

Sq [P) = i^oPo + ^i2^ liPi + mo + — Po + -l^l^Pi + 0{{pQat) , (pia,) ), (34) 

i i 

which yields the following expression for the quark propagator: 

S,{p) = 1 -^loPo--^J:.l.P. + mo + ^-fip^, + vJ:.pl) , ^^^^ 

We consider the tree-level on-shell improvement for this quark propagator requiring that 
Sq{p) should reproduce the form 

Sq{p) = ^ °2°, Y- 2^ ' 2 — - + ('^o pole terms) + 0{al^) (36) 

^q Po + l^iPi + ^i? 

with the appropriate choice for Zq, Zm-, v and ?7, where Zq and Z^, denote renormalization 
factors for the quark wave function and mass defined by 



qR = Z^'q, (37) 

niR = Zmino, (38) 

where mn is the pole mass. It should be reminded that r is a free parameter. We remark 
that terms without the pole in the quark propagator of eq. (p6D yield only contact terms in 
the configuration space, which do not contribute to on-shell quantities in Green functions. 
In terms of the inverse quark propagator, the condition eq. (|36|) is equivalent to 

S^\p) = [Z^ - 2Cmoat] (^^ + mn) + Cat^p" + ml) + 0{at{i]l> + m^)^) + 0{al) (39) 

with C constant. Therefore "(no pole terms)" in eq. ( p6D are not necessary to be 0{a^). 
Comparing the expressions of eqs. (^) and (|36|) , we find at the tree level 

Z-i/2 = 1 - ^moot, (40) 

T 

Z^ = l- -moat, (41) 

z/ = 1, (42) 

r] = l. (43) 

Up to now three types of quark actions with different choices for r and rj have been 
proposed: (i) r = 1 and r] = 1 |^, (ii) r = ^ and t] = 1 |TD[ and (iii) r = 1 and ^ = ^ 0. 



Although the action in the case of (iii) has been most extensively studied numerically, the 
choice of the parameter rj does not meet the condition of eq. ( ^31) that is required from the 
on-shell improvement at the tree-level. This primitive failure makes us consider that it is 
not worthwhile to work on the case (iii) in this paper. We focus on only cases (i) and (ii) 
hereafter. 

Let us first derive the relation between the bare quark mass m,^ and the pole mass m,p. 
Putting Pi = and Pq = imp into the inverse free quark propagator of eq. (p^, the on-shell 
condition yields 



mpQt = log 



m^ottt + r + J{m,QatY + '^rm^^at + 1 
l + r 



(44) 



While in the case (i) with r = 1 we can expand rripat in powers of mo under the condition 
rriQat <^ 1, in the case (ii) with r = C, the condition C,iTLQat = iTiQag <^ 1 is necessary. To 
avoid any confusions we assume moOs ^ 1 from now on. We remark that this assumption 
does not affect any conclusions in this section. 

On the anisotropic lattice we have to be careful about contributions of space doublers. 
Pole masses of space doublers are written as 



m^at = log 
with 



rripat + r + J {raQatY + SrmQaj + 1 



1 + r 



2rN, 



(45) 



< = ^o + i — , (46) 

where Nd components of spatial momentum pi are equal to tt/cs at the edge of the Brillouin 
zone. Although doubler pole masses are always heavier than the physical one irrespective 
of the value of mo, r and ^, their differences in the large limit of ^ are given by 



[ml-mp)as^-Nd-- ^ O [ — ] case (i), (47) 



i l + moflt V4^ 

{rrip - mp)as -^ y^l + 2moa^ + ANd - v^rT2m^ + '^ ( 7 ) ^^^^ *^"''' ^'^^^ 

For the case (i) we find the gap {nip — mp)as diminishes as S, becomes larger. This brings 
a practical problem in numerical studies of heavy quarks: contributions of doublers could 
contaminate signals of hadron states. On the other hand, we are free from this problem in 
the case (ii). 

D. One-loop quark self energy 

The inverse full quark propagator is written as 

Sq^{p) = hoPo + J^i^7iPi + mo + -pI + r]-J2pl - ^P,mo), (49) 



where we take a = rat{;a = at for (i) while a = as for (ii)). One-loop contributions to the 
quark self-energy S(j9, mo) consist of two types of diagrams depicted in Figs. |1] (a) and (b), 
which are expressed by 



T.a{p,m) 



+yt{v.V^ k)S,{p + k)V,%p + k,p) 
+Vi%P,P + k)S,{p + k)V,%p + k,p) 
+V;i(p,P + k)S,{p + k)V,%p + k,p)} GXik) 



(50) 



and 



S.(p) 



d^k 



Y.Y.y2it{p.p)Gt^{k) 



(51) 



A M 



with 



< A;o < -, 

at at 

TT , TT 

<ki<—. 

a. a.. 



(52) 
(53) 



Expanding S(p, mo) = Sa(]9, mo) + S{,(p) in terms of p and mo, we obtain 



S(p,mo) 



167r2 



:C. 



- + ^7oPo(-i^ + S*) + « ^ -iiPi{-L + St) + mo(-4L + S2) 



+ap^((l - 3csw)V2 + <y\) + a^Pi ((1 " 3csw)V2 + ^D 

i 

+amoi7oPo((5 + 3csw)^/2 + (j\) + amo^ Y^ liPi{{^ + 3csw)^/2 + 0-2) 

i 

+am^((5 - 3csw)^ + f^s)] + 0{a]^) (54) 



with L = — log(A^a^) the contribution of the infrared divergence. Here we take cb = ce = 
csw for the clover coefficient. So, S^^''', S2 and ai'^, (J2'*, a^ are independent of a but functions 
of ^, 7] and r parameters. We evaluate these quantities numerically using the Monte Carlo 



integration routine BASES [ITT 



From eqs. (BUp and (El) we obtain 



10 



z-^i^ 



1 + 



9^ C, 
167r2 2 

9" 



(-L + Af) 1 + 



\ 2 167r2 "^ *? , 



l + Tf^C^(3L + A(^) 1 + 



IGtt^ 



""" -i+if^^"^" 



z/ 



1 - iCj^-^i"' - if^^-"»^^"- 

1 _ 2^CfA10) 



167r2 



r? ) 



(55) 

(56) 

(57) 
(58) 



where 



9" ^ So 



167r2 



a 



(59) 
(60) 
(61) 
(62) 
(63) 
(64) 
(65) 
(66) 

We find for tlie anisotropic case tliat the g^aXoga terms disappear in Zq, Z^ for csw =1 as 
well as V and t] for an arbitrary values of Csw, as observed in Zq and Zm for the isotropic 



A(o) 


- y* 




A« 


( t ^2 


^, S2 3(1- Csw). \ 
"^1+ 2 ' 4 ^J' 
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= r (Ti + cTo - 
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case 



12|. Thus Csw = 1 gives the tree level estimate for cb and c^, and we take this value 



for the latter numerical calculation in this section. With this choice for Z^, Z^, v and r/ the 
quark propagator is given by 



Sq{v) 



^-^ 



Po + J2iPi+m 



R 



-hoPo - i Yl ^iPi + ^R 



+«(po + EPi + 



m 



R) 



167r2 V 2 ^2 



Z-'Sf{p) + Z-'a - 



9' 



:C, 



2 ^ 2 , 



(67) 
(68) 



^2 167r2 

where S^{p) is the renormalized quark propagator. We again remark that the term with- 
out the pole in the quark propagator does not contribute to on-shell quantities in Green 



11 



functions. 

We show ^ dependences of one-loop coefficients for Zq, Z^-, y and r] in Figs. 0, |^, | and 
^, respectively. Here we consider both r = 1 and r = ^ cases with r^ = 1, which satisfy 
the tree-level on-shell condition in eq. (0). Although our main concern is whether mog 
corrections could revive at the one- loop level for the r = 1 case, we also present the results 
of the r = ^ case for comparison. Results for A^°\ l^^q^ ■, A^^ and A^^ in the isotropic case 
are already given in Refs. |ll3| , [l^ . They show an agreement with our results with the choice 



of ^ = 1. For the r = 1 case, Figs. ||(b), ^(b) and |^ show approximate linear dependences 
on ^ for A,[^\ A[,^^ and A,|^'^\ which tells us that 0{g^a) contributions to Zm, v and r/ are 
effectively of order g^mas = g'^matC,- We observe similar linear dependences for the r = ^ 
case in Figs. 0(b), |^(b) and |^. From these observations we conclude that mag corrections 
are allowed to revive at the one-loop level. 

This is a reasonable conclusion in view of the on-shell improvement . As far as we know 
there is no symmetry on the anisotropic lattice to prohibit the higher dimensional operators 
with the form of {mas)"'04^^k {n, k > 1), where ^4+^ denotes 4 + k dimensional operators. 
Unless such symmetry is uncovered, the theoretical advantage of the anisotropic lattice over 
the isotropic one would be never confirmed. 

III. ISOTROPIC LATTICE 

In this section we propose a new relativistic approach to control mqa corrections for 
the heavy quarks on the lattice. The idea is based on the on-shell improvement programme 
applied to the heavy quarks on the isotropic lattice. This method allows us to obtain the 
physical quantities in the continuum limit without requiring harsh condition niQa <^ 1 that 
is not achievable in near future. 

Let us consider the general cutoff effects for the heavy quarks on the lattice. Here we 
assume that the heavy quark mass rriQ is much heavier than Aqcd, while the light quark 
mass rrig is lighter than Aqcd- Under this condition we assume that the leading cutoff effects 
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are 



foirriQa) > /i(mQa)aAQCD > f2{mQa) [aAqci)) > ■■, (69) 

where fi{mQa) {i > 0) are smooth and continuous all over the range of rriQa and have Taylor 
expansions at mqa = with sufficiently large convergence radii beyond rriQa = 1. To control 
the scaling violation effects we want to remove the cutoff effects up to /i(mQa)aAQCD by 
adding the counter terms to the lattice quark action with the on-shell improvement. If rriQa 
is small enough, the remaining /2(?7iQa)(aAQCD)^ contributions can be removed by extrap- 
olating the numerical data at several lattice spacings to the continuum limit. Otherwise, 
in case of sufficiently small lattice spacing, the ©((aAqcD)^) errors can be neglected. Our 
aim in this section is to search for the relevant counter terms required in the on-shell im- 
provement and propose a nonperturbative method to determine their coefficients. We also 
show that the behavior of fi{mQa) {i > 1) in the large mqa region can be discussed by 
investigating the static limit. 

A. On-shell improvement on the isotropic lattice 

We first list the allowed operators under the requirement of the gauge, axis permutation 
and other various discrete symmetries on the lattice, where the chiral symmetry is not 



imposed. According to the work of Ref. [|15|, all the operators with dimension up to six are 
given by 

dim.3 :C3(x) = q{x)q(x), (70) 

dim.4 :04(x) = q{x)lpq{x), (71) 

dim.5 -.OU^) = qix)Dlq{x), (72) 

05b{x) = iq{x)o^yFf,^q{x), (73) 

dim.6 :06a(x) = q{x)^^Dlq{x), (74) 

Og,{x) = q{x)Dlpq{x), (75) 

13 



ae(x) = q{x)pDlq{x), (76) 

OU^) = iq{x)j^[D,, F^,]q{x), (77) 

Oeeix) = qix)p'qix), (78) 

C>6fix) = q{x)Tq{x)q{x)Tq{x), (79) 

where F = 1,75,7^,7^75,0"^,^. The higher dimensional operators are related to the lower 
dimensional operators with additional (mga)" corrections with and without classical field 
equations, or otherwise they have the contributions of order (aAqcd)^ or less as the cutoff 
effects. 

These operators lead to a following generic form of the quark action on the isotropic 
lattice: 



s'r = E 



C303{X) + C^O^ix) + Y. C5^05^{X) + J^ Ce^Oe^{x) 

i=a,b i=a,...,f 



(80) 



where C3, . . . ,C6/ are functions of the bare gauge coupling g and the power corrections of 
ma. We first remark that the ma corrections to the quark mass term and the kinetic term 
can be absorbed in the renormalizations of the quark mass Zm and the wave function Zq. 
For the sake of convenience we choose C3 = m,Q and C4 = 1. 

In the next step we reduce the number of basis operators with the aid of the classical 
field equations. It is easily found that O^a, C^Qb, Oqc and Oge can be related to the quark 
mass term or the kinetic term. In the on-shell improvement these operators are redundant 



and can be eliminated from the action of eq. (|80D . The operator O^a, however, is used to 
avoid the species doubling and the value of its coefficient c^a is given by hand. 

The remaining operators are O^b, C'ea, C'ed and Oe/, whose contributions as the cutoff 
effects are estimated in Table |I[ We use the classical field equation for g(x)7o-Dog(x) in Oga- 
Here it should be noted that Aqcd means the order Aqcd or less throughout this paper. In 
some cases the actual contribution may become smaller. For example, the contributions of 
the bilinear terms whose Dirac matrices consist of off-diagonal components are suppressed 
by the extra Aqcd/'^q for the heavy quarks compared to the light quarks. The operator 

14 



O^b is the so-called clover term, for which the nonperturbative method to determine the 
coefficient c^t, in the massless limit is already established ||T6|. However, the contributions 
of {ma)^05b {n > 1) cannot be neglected in the present condition that allows rriQa ~ 0(1). 
For O(aAQCD) improvement the coefficient c^b has to be adjusted in the mass dependent 
way. The differences in magnitude between the time and space components in O^a originate 
from the violation of rotational symmetry on the lattice with finite lattice spacing. While 
the contributions of the space components are found to be negligible, those of the time 
components should be removed. We also find the contributions of the four-quark operators 
in O^f are negligible. 

The generalization of the above argument to any operators with higher dimensions makes 
the discussion more transparent. Let us consider an arbitrary operator with A + k dimension, 
a^04^^k, where we write the lattice spacing a explicitly. The operator (^4+^ contains / pairs 
of q and q and n covariant derivatives D^ with 4 + k = 3 x I + n. Using the classical 
field equation, some (but not all) of covariant derivatives can be replaced by the quark 
mass m. For / > 2 the largest possible power of the scaling violation is (?7ia)"(aAQCD)^'~^- 
Therefore the operators which contain four or more quarks are irrelevant for the ©(qAqcd) 
improvement. All the relevant contributions come from the quark bilinear operators. With 
the aid of the classical field equations, they can be reduced to 

{maYa~^q{x)q{x) (81) 

(ma)"-^g(x)7o/^og(x), (ma)"-^ ^g(x)7iAg(x) (82) 

i 

{maY-'^aq{x)Dlq{x), {maY-^aY,Q{x)Dlq{x) (83) 



(ma)^ "^ai 



^q{x)aQiFoiq{x), {ma)"- '^ai^q{x)aijFijq{x), (84) 



for n > 0. The time and space components of O4 and O^a^b should be treated separately in 
case of finite ma, where the space-time asymmetry reflects the contributions of the higher 
dimensional operators that break the rotational symmetry. Now we know that the seven 
operators are needed for the ©(oAqcd) improvement. Since three coefficients among these 
seven operators can be absorbed in Z^, Zq and the Wilson parameter r^ for the time deriva- 
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tive as already explained, the remaining four coefficients have to be actually tuned. 
In conclusion, at all order of ma, the generic quark action is written as 

rta_, 



riimp \ ^ 



mQq{x)q{x) + q{x)-iQDQq{x) + v^q{x)-iiDiq{x) - —q{x)Dlq{x) 



i i i,j 



(85) 

where we are allowed to choose r^ = 1 and the four parameters z/, r^, ce and cb are to 
be adjusted. In general these parameters have the form that X = J2n^n{g'^){'m'0')"' with 
X = u, Ts, Ce and cb, and Xq should agree with the one in the massless 0{a) improved 
theory: uq = 1, {rs)o = ^t = 1, {ce)q = {cb)q = Csw 0]- Note that v = 1 + 0{{maY) and 
rs = Tt + 0{ma) since the space-time asymmetry arises from Lorentz non-covariant terms 
such as O^a via the on-shell reduction, accompanied by extra {maY factors. 

From the above consideration, the leading scaling violation in the massive theory, except 
for E^=iQ''"((7', log a) (mQa)" in Z^ and Z^, is E^=oCf (^', loga)(mQa)"aAQCD for the 
Wilson quark action, or I]J^i C^^{g^, loga)(mQa)"'aAQCD for the (massless) 0{a) improved 
SW quark action. Here we should notice that the contribution of (i^—l) J2i (i{.x)liDiq{x) is of 
order roQa^ Kqqj) / mq ~ mqa? Kqqd ~ aAqco in the heavy quark region. This implies that 
once we fix the pole mass from some spectral quantity, the cutoff effects in other spectral 
quantities are at most of order oAqcd, not {mQaY, for the Wilson and the (massless) 0{a) 
improved SW quark actions. It should be noted that the quark wave function does not affect 
on the spectral quantities. If z/, r^, ce and cb are properly adjusted in the mass dependent 
way, the remaining scaling violations are reduced to Y,'^=QC^'^^{g'^ ,\oga){mQaY{aKq^cTi)'^ ■ 

This relativistic argument about the on-shell improvement on the massive quarks helps 
us understand some numerical results for heavy quark physics previously obtained by using 
the Wilson and SW quark actions under the condition mqa ~ 0(1). We find a good 
example in Fig. 2 of Ref. |T^, which compares the difference between twice of the heavy- 
light meson mass and the heavy- heavy one, 2mHL — rnnH, obtained by using the pole mass 
or the kinetic mass defined by d'^Euu/dpf. We observe that 2mBL — i^hh with the pole 
mass is consistent with the experimental values within the ambiguities of order oAqcd as 
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expected, while 2mHL — i^hh with the kinetic mass gets deviated from the experimental 
values further and further as niQa becomes larger. Recall that the authors of Ref. [Q suggest 
from a nonrelativistic point of view that the kinetic mass should be used for analyses on the 
heavy quark quantities. Our relativistic argument, however, tells us that the use of the pole 
mass makes the remaining cutoff effects O(aAQCD)- Since the difference between the kinetic 
mass and the pole mass is of order (mga)^ at the tree- level, the use of the kinetic mass 
eventually yields unwanted additional (rriga)^ errors. This is the reason why the results 
of 2mHL — ^HH with the kinetic mass show considerable deviation from the experimental 
values. 

At the end of this subsection we have to remark one point. The ma corrections in 
Zjn and Zq, though they are irrelevant for spectral quantities, become important, together 
with other renormalization factors in case of calculating the quark masses and the matrix 
elements of various composite operators. In Sec. [1IIF| we will show how to calculate these 



renormalization factors nonperturbatively including the {rriQa)"' corrections. 

B. Improvement of the axial current 

The cutoff effects in the correlation functions of local composite fields are originated 
from not only the action but also the composite fields themselves. In this subsection we 
demonstrate the on-shell improvement on the axial current, which is relevant for calculation 
of the heavy-light pseudoscalar meson decay constants like /b and fo- 

The axial current in the Nf flavor space is given by 

A;{x)=qix)^^^,-qix), (86) 

where A"/2 (a = 1, . . . , A^? — 1) are generators of SU (Nf). The improvement of this op- 
erator is performed in the same way as the improvement of the quark action. After some 
consideration we find that it is sufficient to consider only the dimension four operators for 
the ©(oAqcd) improvement: the higher dimensional operators can be reduced to ^JJ or the 
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dimension four operators multiplied by (ma)"' using the classical field equations, or other- 
wise their contributions are of order (oAqcd)^- The requirement of various symmetries on 
the lattice allows the following dimension four operators: 

dim.4 iAAaT^ix) = q{x)i,,D^—q{x) + g(x)S^75yg(x), (87) 

{Aib)l{x) = q{x)-i^a^^D^—q{x) - q{x)D^a^y-i^—q{x), (88) 

where we do not take the sum on the index /i. We find that (^4?,)^ is related to {A^aY^ 
and A'^^ with the aid of the classical field equations, which means the operator {A^hY^ is 
redundant for the ©(oAqcd) improvement. 

After all the improved axial current is written as 

{Air^{x) = A^ix) + d,{g\ma){A,a)l{^). (§9) 

The parameter (i^ has been already calculated in the massless case. Perturbative estimate 
gives df^{g'^,ma = 0) = —0.00756(7^ [|l^, which is fairly small in magnitude. From this fact 
we think that it would be sufficient to evaluate c?^ perturbatively in the massive case. The 
lattice operator (.4.")'™'' is related to the continuum operator with the renormalization factor 
Za- 

(^«)-°(x) = Z^{g',\oga,ma){A';y^^{x). (90) 

The power corrections of ma in Za need to be under control to obtain the heavy-light 
pseudoscalar meson decay constants defined in the continuum regularization scheme. In 
Sec. [1111^] we will show how to remove the ma corrections in Za with a nonperturbative 
method. 

C. Benefit of chiral symmetry 

Although the above discussions are free from the chiral symmetry, it is also interesting to 
look into what can be changed by the presence of the chiral symmetry. For the convenience 
we treat the quark mass matrix in the Nf fiavor space 



M 



mi 

1712 










(91) 



y u u ■ ■ ■ niN; J 
as a spurious field M. which transforms hke 

M ^ VrMVI 



(92) 
(93) 



under SU(A^/)lxSU(A^/)/j, where Vl and Vr are elements of the fundamental representation 
of SU(A''j)i and SU(A^/)i? respectively. In terms of quark fields q and g, SU(A^/)l and 
SU(A^/)ij act on the left and right handed components, 



1-75 

-I + 75 

qi = ?^— ' 

1+75 
QR = 2 Q: 

_1 - 75 
Qr = Q — 7^ — , 



whose transformation properties are given by 



(94) 
(95) 
(96) 
(97) 



qi- 


-^ VigL, 


Ql- 


- QlVI 


Qr - 


-* VRqR, 


Qr - 


-^ QrVr- 



(98) 
(99) 

(100) 
(101) 

We work all the calculations assuming this symmetry and at the end of calculations we can 
choose M = M"^ = M. 

With this artifice let us consider which operators among O3, . . . , Oqj in the quark action 
are allowed under SU(A'^y)LxSU(A^/)ij symmetry. We easily find that the dimension three 
and five operators are not allowed. As for the dimension six operators, some four-fermi 
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operators are excluded. We also observe that the power corrections of mqa emerge as the 
form 

M=^" ■ (SU(A^/)l X SU(A^/)ij invariant operators) (102) 

with n > 0, which means 



') 



g(a;)M2"+ig(x), (103) 

q{x)M^^^,D^q{x), (104) 

iq{x)M^^+'a^,F^,q{x) (105) 



are allowed, while 



q{x)M^''q{x), (106) 

q{x)M'''+'-f,D^q{x), (107) 

2g(x)M=^V^,F^,g(x) (108) 

are forbidden. This may be advantageous in controlling the cutoff effects as mqa becomes 
smaller away from one. Even for the chiral non-invariant quark action such as the SW 
quark action, however, the leading cutoff effects except for Zg and Zm are (mQa)"aAQCD 
with n ^ 0, which are of the same order as those in the chirally symmetric actions, once the 
coefficient of O^f, in the quark action is nonperturbatively tuned in the massless limit. 

As for the improvement of the axial current, the similar argument can be applied. The 
dimension four operators {A^ajfi and (^4^)^ are not allowed by the chiral symmetry and the 
power corrections of rriQa are restricted to the form of eq. ( |102|) . As an example, 

qix)M'-+'^,D^,^qix) + qix)M^^^+''D,^,^qix), (109) 

with n > are allowed. 
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D. mqa corrections at tree-level 

In order to control the rriQa corrections it should be essential to nonperturbatively de- 
termine the renormalization factors Z^ and Zq and the four parameters z/, r^, ce and cb- 
However, we think it is instructive to first investigate the mqa corrections at the tree-level. 

Zq, Zm, i^ and r^ can be determined by demanding that the tree-level quark propagator 
Sq{p) derived from eq. (|85| ) should reproduce the relativistic form 



SqiPo,Pi) = -^ °2°i Y^ 2 I ' 2 — ~ + ("^^ P°^^ terms) +0((pia)^) (110) 



Zq Po + EiPf + mt 



around the pole. Imposing pi = we first obtain 



rrip = log 


mo + rt + ^ml + ^nnio + 1 

i + n 


mo 


Zq = cosh(mp) -I- r(Sinh(mp). 


We then find with finite spatial momenta 


sinh(mp) 


nip 


cosh 
rs = 


'nip) + risinh(mp) sinh(mp) 
nip nil 



(111) 

(112) 
(113) 



(114) 

(115) 

We should notice that the niqa corrections start at 0{{niqaY), not 0{niqa), in the u 
parameter as expected. Figure ^illustrates the nipa dependences of Zq, Z^, y and r^ in case 
oi Tt = 1. We observe that the nipa dependences of u is relatively mild compared to those 
of Zq, Zm and r^. 

To fix the ce and cb parameters we consider the quark-quark scattering amplitude de- 
picted in Fig. 1^. The improvement condition is that ce and cb should be chosen to reproduce 
the following form of scattering amplitude at the on-shell point removing the niqa correc- 
tions, 

T = -g^u{p')-ff,u{p)Df,^{p - p')u{q')-i^u{q) 

-g\{q')-f^u{p)D,,{p - q'Mp'hMq) + ^((p.a)^, (g,a)^ {p',a)\ {q'.af), (116) 
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where Dfj^^, is the gluon propagator on the lattice. Notice that with the use of the Gordon 
identity (; on-shell condition for external spinors u, u) the quark-gluon interaction induced 
by the clover term can be transformed into the ordinary quark-gluon vertex: 



1 



= u{p')-[ism{p'Q) + ism{po) + 7o{2(mo + n) - rt^cos^p'o) + cos(po)) 

+r, 5^(2 - cosip'i) - cosipi))}]uip), (117) 

I 

u{p')iaio{sm{pQ) - sm{po))u{p) + u{p') ^ i'iaii{?,in{p[) - sm{pi))u{p) 

I 

= u{p')[u{ism{p[) + ism{pi)) + 7i{2(mo + n) - rt{cos{po) + cos(po)) 

+r, ^(2 - cos(p;) - cos{pi))}]u{p). (118) 

I 

This improvement procedure with the finite quark mass is an extension of the previous work 

||T9| that determined the csw = ce = cb parameter up to one-loop level in the massless case. 

After some algebra with the aid of eqs. (|117|) and ( p.l8|) , we obtain 



Ce = niy, (119) 

Cb = Ts, (120) 

where z/ and r^ are already determined from the on-shell improvement on the quark propa- 
gator. 

You may have already noticed that our values for u, r^ and ce are different from those 
derived in Ref. P|. This difference originates from whether the on-shell improvement is 
implemented in the relativistic way or in the nonrelativistic way. (more precisely, in case 
that the Lagrangian does not retain the rotational invariance on the Euclidean space-time, we 
need both the momentum operator and the Hamiltonian to discuss the rotational symmetry.) 
For example, both methods give the same relation for the v and r^ parameters from the 
dispersion relation: 

2 sinh(mp) 

V + rssmh(mp) = {cosh(mp) + rtsinh(mp)} . (121) 

trip 
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In the nonrelativistic approach, however, v and r^ are not distinguishable due to the lack of 
relativistic informations. Generally speaking, we do not have sufficient number of nonrela- 
tivistic conditions to fix the coefficients of the relativistic operators with higher dimensions 
because the degrees of freedom of quarks are smaller in the nonrelativistic approximation 
compared to the relativistic case. Although the nonrelativistic approach with the Wilson 
type quark action ^ has been considered to work better than the NRQCD in the charm 
quark region where the sizable relativistic effects are expected, this is not necessarily assured 
because this approach does not meet the relativistic on-shell improvement. 

E. Large ttiqo and static limit 



Although we have restricted ourselves to the case of finite mqa so far, it is worthwhile 
to show that we can derive the static quark action from eq. (pSf ) by taking mqa ^ oo. In 
terms of the heavy quark field h[x) defined by 

AEt 



q[x) 



h{x), 



(122) 



where E = ±imp, the lattice action in eq. (|85D becomes 

cimp \ ^ 



ft 3rs \ ^ , , . , , e*^ 



1 + ^ + ^^ ) h{x)h{x) - - — h{x){rt - 7o)f/o(x)/i(x + 0) 
mo rriQ J Ztuq 



-iE 



2m, 



-h{x + 6){rt + -fo)Uo{xyh{x) 



+ 



[iy,CE,rs,CBy 
mo 



Taking m,o -^ oo and setting r^ = 1 for simplicity, we obtain 

^imp ^ ^ \h{x)h{x) - h{x + 6)^^^^Uo{xyh{x) 
a; L -^ 

for the heavy quark(ii^ = im,p), or 

E 

X 

for the heavy anti-quark(£' = —imp), since 



Qimp ^ 



h{x)h{x) - h{x)^-—^Uo{x)h{x + 0) 



(123) 



(124) 



(125) 
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Jl = i£.0^ , (126) 

mo mo yrripj 

^ = ^-o(^], (127) 

mo mo ympj 



Q-irip 



r-^ 



1, (128) 



mo 
where m.p = log|l + m.o\ is used. We replace ^-^h{x) = h{x) from the property that 
7o/i(a;) = h{x) for the quark and '~foh{x) = —h{x) for the anti-quark. Thus we exactly obtain 
the static quark action, where the quark moves forward, or the static anti-quark action, 
where the anti-quark moves backward in time. 

The existence of the static limit may give some constraints on the mass dependent scaling 
violation for m,Qa ^ 1. As an explicit example, we take the heavy-light pseudoscalar meson 
decay constant J'hl- One can extract J'hl from the correlation function of the heavy-light 
current Ao^{x) = qH{x)'~^olb(lL{x) with zero spatial momentum following 

Zf\A^\t)A^'^m ^ flL^HLe~^"'\ (129) 

where Z'^^ denotes the renormalization factor of the quark wave function at the tree- 
level. On the other hand, in the static limit, the current is given by h{x)'~fol5(lL{x) = 
^yrn^e"^'^^ Aq ^ (x) , and the correlation function 

mo(^?^(t)^^^(0)) ~ C'e-^^* (130) 

has a well-defined limit. From the relation ZJ^^^ = e™'^ ^ mo with r^ = 1 and m.^^ — 
ntQ + Ai?, we obtain 

fHL = -^. (131) 

y/rriHL 

Since the continuum limit can be taken in the static theory, C should behave as 

^ Wo + WkiaAQCD)" + 0((oAqcd)'+') (132) 



A 3/2 

where A; = 1 for the Wilson light quark action or A; = 2 for the (nonperturbatively tuned) 
SW light quark action. Therefore the following relation holds for m,Qa -^ oo: 
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Jhl _ / Aqcd 

AqcD V ^HL 



Wo + Wk{aAQCDr + 0{{aAciCBr') . (133) 



On the other hand, the decay constant should behave 



Aqcd 



^o(AqcdMq) 

1 + v^iaAQCD)" + v^{mQa){aAciCBr + 0((aAQCD)'+\ (aAQCD)"+')l (134) 



X 



from our consideration in the previous sections, where vj^ is the scahng violation caused by 
the light quark action, and v^{mQa) conies mainly from the Wilson/SW heavy quark action 
for n = 1 or from the action of eq. (^) with the nonperturbatively tuned u, Vg, ce and cb 
for n = 2. Comparing eqs. (|133|) and ( |134|) , we find 



rriHL 
and 



Vo{Aqci) /mq) -> J-^^Wo (135) 



v^- 


-^Wk, 


Vnimqa) - 


-^0 


ioT n < k 


(136) 


Vk+^ni^Qa)- 


^Wk, 


for n = A; 






(137) 


vl- 


■^Wk, 


V^i^QO') - 


^Wn 


ior n > k 


(138) 



in the limit rriQa -^ oo, where Wk and Wn are universal constants independent of the choice 
of the heavy quark action. Note that we expect Vq can be expanded in terms of of mq/ Aqcd 
near the chiral limit. 

The important point is that the function v^{mQa) becomes constant (or even vanishes) 
in the limit mqa -^ oo. On the other hand, we also know that the behavior of v^{mQa) 
is benign under the condition mQa<^l. From these facts it would be reasonable to assume 
that the function v^{mQa) behaves modestly all over the range of mqa from the massless 
limit to the static limit. This assumption is supported by previous numerical results for the 
heavy-hght decay constants fo^ and fs^ which are obtained by using the pole mass under 
the condition mqa'^l with the Wilson quark action. Figures 4 and 5 in Ref. |T^ show that 
the naive continuum extrapolation with a linear form for fa^ and fs^ gives "reasonable" 
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values in the continuum limit. This outcome is understood as follows: Since v^{mQa) is 
a modest function in terms of mqa, the leading scaling violation effects v^ {mQa)aAQCD 
for fos and Jb^ are effectively removed by the naive linear extrapolation. However, the 
O(aAQCD) cutoff effects cannot be completely removed. The remnant is estimated to be 
\vi i^TT-Qdmax) — v^ {mQajam)\(^iw-^QC'D, whcrc Omax and a^am a^c the maximum and minimum 
of the lattice spacing used for the continuum extrapolation, and Oav = amaxflmin/lcmax— Qmin)- 
Although we have focused on the case of heavy-light pseudoscalar meson decay constants, 
the above arguments on scaling violation effects can be easily generalized to any observables 
which can be defined in the static limit. 

F. Nonperturbative renormalization 

Let us turn to a nonperturbative determination of Zg, Z^, v-, Ts, cb and cg. We first 
consider the v and r^ parameters. Since the rotational symmetry breaking due to the mqa 
corrections deviates the v and r^ parameters from one, it would be a reasonable way to adjust 
them such that the correct dispersion relations are reproduced for some hadronic states. We 
think a set of dispersion relations for the heavy-heavy and heavy-light mesons is a good 
choice. A previous study demonstrated a clear distinction between the the kinetic masses 
defined by d'^E-nn/dpf for the heavy-heavy meson and d'^Em^/dpf for the heavy-light meson 
(see Fig. 1 in Ref. |l^), which tells us that the dispersion relations for the heavy- heavy and 
heavy-light mesons give two independent conditions to fix both u and r^ parameters. We 
point out that to avoid the ambiguities coming from the clover term it would be better to 
consider the dispersion relations of the spin averaged meson states over the pseudoscalar 
and vector channels. This is motivated by an observation that the clover term causes the 
hyperfine splitting. 

Nonperturbative determination of cb and ce is a little bit troublesome. Although in the 
massless case the clover coefficient can be determined with the aid of the PCAC relation, this 
method does not work with the massive case. The reason is that the chiral symmetry allows 
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the clover terms with the odd power of ma corrections {maY^^^iq{x)(j^pFf^yq{x) {n > 0) as 
discussed in Sec. [111C|. This imphes that even the chirally symmetric quark actions, e.g., 



the domain wall and the overlap quark actions, suffer from this difficulty. However, we can 
at least evaluate cb and ce perturbatively up to one-loop level by extending the calculation 
in Sec. [111D| . In this case the remaining cutoff effects are of order a^aAqcD, which might be 
small enough for numerical studies. 

As for a nonperturbative determination of renormalization factors, we consider the use 



of the Schrodinger functional (SF) method |^0||. The renormalization of the quark mass is 



made through the renormalizations of the axial current and the pseudoscalar density using 
the PCAC relation: 

m( = l/L) = ^^(^' ^^^ 9, ma)df,q{x)-ff,-f5q{x) ,^^^. 

Zp{g,L,ThL,g,ma)q{x)'y5q{x) 

where L is the physical box size and g and rh are the renormalized coupling and quark mass 
in the SF scheme. Za and Zp with the finite quark mass are defined by 

ZA{g, rfiL, g, ma) = —- — "" ]. , (140) 

fA{xo = L/2) 

Zp{g, L, ml, g, ma) = -— " \ , (141) 

Jp\Xq — Lj l) 

where c^ is a mass-dependent constant, and /^, fp and /i are the correlation functions 
given by 

/-4(xo) = -\j d?0z{q{x)^,^,q{x)C{y)l,az)), (142) 

fp{xo) = -IJ d^yd'mx)i5q{x)ay)i5m), (143) 

1 



/i = -^ J d'ud'vd'yd'z{C'iu)^,C'iv)ay)l5Cm (144) 

with (, (', ( and (' the boundary quark fields. We illustrate the function fA,p and /i in 
Fig. Ij. As for the boundary conditions we refer to the description in Ref. |^ . 

Although each of Za, Zp, (^(x)7^75g(x)) and {q{x)'~^^q{x)) has the power corrections of 
m,a, they are canceled out in the combinations Z^(g(x)7^75g(x)) and Zp{q{x)'^^q{x)) and 
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the 0((aAQCD)^) uncertainties are left. This assures us to take the continuum hmit for 
the renormahzed matrix elements and also for the renormalized quark mass of eq. ( |139D . 
These quantities, however, are defined in the SF scheme with the finite quark mass and 
hence different from those renormalized in the SF scheme at the massless point, which are 
exactly what we want. Therefore, we still need the finite renormalization factor to make the 
conversion between the two renormalization descriptions. This can be obtained by taking 
the continuum limit of the ratio 

Za^Q, friL = 0, g', ma' = 0) 



ZAig,mL,g',ma') 
Zp{g, L, rfiL = 0, g', ma' = 0) 



(145) 
(146) 



Zp{g,L,mL,g',ma') 

for several g' chosen to satisfy m,a' <^ 1. We remark that the physical size L for the SF 
scheme can be taken to be much smaller than that for the measurement of the spectral 
quantities and the various matrix elements, allowing us to access the finer lattice spacing 
with g'. 

IV. CONCLUSIONS 

In this paper we have first examined the validity of the idea of anisotropic lattice by 
making a perturbative calculation of the quark self energy up to 0{g^a). Our results show 
that the m^qag corrections revive through one-loop diagram. We also find that on the 
anisotropic lattice the four parameters must be adjusted to remove all the terms of order 
a even in the massless case. From a theoretical point of view the anisotropic lattice is not 
necessarily advantageous over the isotropic one. 

In the second part of this paper we have presented a new relativistic approach to the 
heavy quarks on the lattice. The idea is based on the relativistic on-shell improvement 
with the finite m^ga corrections. We have shown that the cutoff effects can be reduced to 
©((oAqcd)^) putting the {mqa)'^ corrections on the renormalization factors of the quark 
mass Zm and wave function Zq. As far as the spectral quantities such as hadron masses 
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are concerned, the (mQa)"' corrections in Zm and Zg do not matter: Zg does not affects on 
the spectral quantities and the {rriQa)"' corrections in Z^ can be handled by employing the 
pole mass fixed from some spectral quantity. On the other hand, in case of calculating the 
quark mass or the various hadron matrix elements we can control them by determining the 
renormalization factors nonperturbatively. 

The leading scaling violation for various types of actions are summarized as follows. 
/i(mQa)aAQCD for the Wilson quark (/i(0) 7^ 0) and the 0{a) improved SW quark 
(/i(0) = 0), while /2(mQa)(aAQCD)^ for our proposed action with mass-dependently tuned 
z/, Ts, ce and cb- Therefore, if the magnitude of fi{mQa) is 0(1), the scaling violation for 
heavy hadron masses with the 0{a) improved SW quark might be as bad as that for light 
hadron masses with the ordinary Wilson quark. For sufficiently small rriQa, we can remove 
the leading scaling violations by extrapolating the data at several lattice spacings to the 
continuum limit. Even if mga ~ 0(1), /2(mQa)(aAQCD)^ in our proposed action is expected 
to be negligibly small in case of (oAqcd)^ ~ 0.01 

Our relativistic approach has the strong point over the nonrelativistic ones: the finer 
the lattice spacing becomes, the better the approach works. This is a desirable feature 
because we can take the full advantage of configurations with finer lattice spacing generated 
to control the cutoff effects on the light hadron physics. We are going to perform a numerical 
simulation to test the ideas presented in this paper. 
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FIG. 1. One- loop diagrams for the quark self energy. 
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FIG. 2. ^ dependences of (a) Aq and (b) Aq in the renormalization constant of the quark 
wave function, ce and cb are chosen to be 1. Errors are within symbols. 
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FIG. 3. £, dependences of (a) Am, and (b) Am in the renornialization constant of the quark 
mass, ce and cb are chosen to be 1. Errors are within symbols. 
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FIG. 4. ^ dependences of (a) Au and (b) Ajy in the u parameter, ce and cb are chosen to 
be 1. Errors are within symbols. 
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FIG. 5. ^ dependences of Ar; in the r] parameter, ce and cb are chosen to be 1. Errors are 
within symbols. 
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FIG. 6. Tree-level values for Z„, Z^, i^ and r^ as functions of vtipa. We choose rt = 1. 
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FIG. 7. Tree- level diagrams for the quark-quark scattering. 
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FIG. 8. Quark diagrams contributing to (a) fA,p{x) and (b) /i. Oa,p are inserted at the point 
X. C and C denote the boundary conditions for the gauge fields at t = and T. 
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TABLES 
TABLE L Expected magnitude of the cutoff effects due to the higher dimensional operators in 
the quark action. niQ denote the heavy quark masses. 

operator hght heavy 

Csb : iJ2iQi.x)aoiFoiq{x) aAqcD oA^^oMq 

iT,ijQ{x)(^ijFij<l{x) aAqcD oAqcd 

Cea : q{xhoD^q{x) (oAqcd)^ (.mQafq{x)q{x), {mqaf q{x)-fiDiq{x) , {mQa)q{x)Dfq{x) 

q{x)jiDfq{x) (aAgcD)^ (aAgcD)^ 

Oed- iT,iQ{xhi[Do,Fio]q{x) (oAqcd)^ a^A^c^/mQ 

iT,iQ{xho[Di,Foi]q{x) (cAqcd)^ {aAQcuf 

iJ2iQix)-fi[Dj,Fij]q{x) (oAqcd)^ o^A^cdMq 

Oef : (oAqcd)^ («Aqcd)^ 
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